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TURBULENT SOLUTION OF THE NAVIER-STOKES EQUATIONS 
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R. G. Deissler 
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Cleveland, Ohio 44135 

ABSTRACT 

To study the nonlinear physics of uniform turbulent shear flow, the unav- 
eraged Navier-Stokes equations are solved numerically. This extends our pre- 
vious work in which mean gradients were absent. For initial conditions, modi- 
fied three-dimensional-cosine velocity fluctuations are used. The boundary 
conditions are modified periodic conditions on a stationary three-dimensional 
numerical grid. A uniform mean shear is superimposed on the initial and Doun- 
dary conditions. The three components of the mean-square velocity fluctua- 
tions are initially equal for the conditions chosen. As in the case of no 
shear the initially nonranoom flow develops into an apparently random turbu- 
lence at higfier Reynolds numuer. Thus, randomness or turbulence can appar- 
ently arise as a consequence of the structure of tne Navier-Stokes equations. 
Except fur an initial period of adjustment, all fluctuating components grow 
wuh time. The initial equality of the three intensity components is de- 
stroyed by the shear, tne transverse components becoming smaller tnan the 
longitudinal one, in agreement with experiment. Also, the shear creates a 
small-scale structure in the turbulence. The nonlinear solutions are comparea 


witn linearizeo ones. 


INTRODUCTION 


This is an extension of our work in which the development of isotropic 
turbulence was examined numerically.^ In that work the nonlinear transfer 
of energy to smaller scales of motion and turbulent dissipation were studied. 

Another important process is the production of turbulence by a mean 
shear. Most turbulent flows, both those occurring in nature and those which 
are manmade, are in fact shear flows, where the turbulence is produced and 
maintained by the shear. Becaise of the added complexity, the nonlinear prob- 
lem of turbulent shear flow is even more difficult than that of isotropic tur- 
bulence. So it is not surprising that little progress iias been maae in ob- 
taining an analytical solution from first principles. An attempt to obtain a 
numerical solution would soem to ue in order. 

Conceptually, tiie simplest turbulent snear flow (although certainly not 

the simplest to produce experimental ly‘‘ ) is one in which the turbulence is 

homogeneous and uniformly sheared. At least two significant numerical studies 

j 4 

of that type of turbulence iiave recently been made. ’ In both of tnose 
studies random initial conditions witti a range of eddy sizes were used. 

In tlie spirit of Kef. 1, the present numerical study of un if ormly-sheareo 
turbulence starts with simple determinate initial conditions winch possess a 
single length scale. In this way we can study how tne turoulence develops 
T rom nonturbulent initial conditions. Moreover much higher Keyrio I os-numoer 
riows can be calcalateu wild a given numerical grio when a single length scale 
is initially present, at least for early and moderate times. Some results 
using a perturbation senes were outaineo in Kefs. t> ana o, but tne calcula- 
tions could not be carried viiry tar in tiiiie. Moreover tne compotier ts ot the 
initial intensities were not ('‘quai, as ttiev are litre. 


As will be seen, several interesting results which could not be obtained 
in the previous work on turbi'lent shear flow are obtained in the presetit 
study. One of the significant findings is that the structure of the turbu- 
lence produced in tite presence oi a strong shear is much finer than that pro- 
duced in its absence. 

THE INEAk problem 

As in Ref. i, tlie Navier-Stokes anc continuity equations for an incompres- 
sible fluid are written in dimensionless form as 


au a( u u ) ~ d^u . 

i_ ' 1 K ’ _ ap ^ 1 

9t^ 3x, ax- ax, ax, 
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Neie tnai. the stars on eiiriens ionu i quaniities are oiiiitteo for curresponcn ng 
u iinension less quaniities, liie subscripts can take on the values 1, 2, ana 3, 
ano a rep^.w-itea sui)script in a te'rm muicaies a summation, [tie quantity 
Li ^ IS an i ns taut arunmjs velocity component, x. is a space cooroifidte, x,^ is 

'k 

a character ] si 1 k. ienqtii, t is the time, p is the density, v is the 

* 

KinemaliL viscnsiiy, and p is tne instantaneous pressure. In oroer to 
c)hlaifi an t^xpiicii fM]ualu.ni tor Uie pressure, we idKO tne oivergence ot 
[A], [ Ia ana ap^iiv tiit: continuity tq. (C) to get 
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(3) 


In the remainder of the analysis it Mill be found convenient to use the set of 
Eqs. (1) and (3), rather than Eqs. (1) and (2). 

The expression assumed for the initial disturbance is, in dimensionless 
form 


where 



dU 


1 


cos Xr 
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and 


Xq * 
U . = — U.. 

1 V 1 


The quantity a*'’ is an initial velocity amplitude or Fourier coefficient 
of the disturbance, is an initial wave number vector, U^* is a mean 

velocity component, and is the Kronecker delta (equals 1 for i » j 

♦ ♦ >n n 

0 for i ^ j). The quantities q"*x are dot products (q *x * qj x^ 


and 


+ q^ X 2 + q^ x^). Equation (4) reduces to the initial condition in Ref. 1 
for dU^/dx 2 = 0 (no shear). In order to satisfy the continuity condition, 
Eq. (2), 


n^n 

a^di 


0 . 


For the present work, as in Ref. 1, we set 


(5) 


4 


( 6 ) 


- k{2. 1. 1). . k(l, 2. 1). a^ . k (1. 1. 2). 

q\ - (-1, 1, 1). a (1, -1, 1), and - (1, 1, -1), 

where k is a quantity that fixes the Reynolds number. In addition to satis- 
fying continuity, Eqs. (6) give 


T T T ■? 

ui = U2 . U3 - Uq 


( 7 ) 


at the initial time, where u^ s Uj and the overbars indicate 

averaged values. Thus, Eqs. (4) and (6) give a particularly simple initial 

condition, in that we need be concerned with only one component of the mean- 

square velocity. Note that it is necessary to have at least three terms in 

the summation in Eq. (4) to satisfy Eq. (7). 

In order to carry out the numerical solution of Eqs. (1) and (3) subject 

to initial conditions (4) and (6), we use a stationary cubical grid with a 
3 

maximum of 32 points and with faces at x^. = 0 and 2w. For boundary 
conditions we assume modified periodicity. That is, we let 


= '"i'x-b, * *n «J2 '«l''“>'2 


J J 


(8) 


and 


^*^^x -=2ir+b . ' 

0 J «3 

for any bj. Note that Eq. (8) is not a tensor equation. It is the usual 
periodicity condition, but with a superimposed mean shear, and is consistent 
with the initial condition given by Eq. (4). Equation (8) is used to calcu- 
late numerical derivatives at the boundaries. 


5 


The spaclal- and time-differencing schemes are essentially those used by 
Clark et al.^ and In Ref. 1. That Is, for the spaclal derivatives In 

Q 

Eqs. (1) and (3) we use centered fourtiv-order difference expressions. For 
time-differencing we use a predictor-corrector method with a second-order 

Q 

(leap-frog) predictor and a third-order (Adams-Moulton) corrector. The 

Poisson equation for the pressure (Eq. (3)) Is solved directly (no Iteration) 

by a fast Fourler-transform method. This method preserves continuity quite 
i 

well (v*u a 0). 

Some of the results will be extrapolated to zero mesh size in an effort to 
obtain more accuracy. The method of fourth-order extrapolation In Ref. 1, 
which is consistent with fourth-order differencing, is used here. However to 
increase the accuracy, more of the grid-point spacings are chosen close to 
zero: 32^, 24^, and 16^ grid points, rather than the 32^, 16^, and 8^ points 
in Ref. 1 are used. Equation (10) in Ref. 1 must then be replaced by 

t^ a 1.6814 t^ - 0.7182 + 0.0368 t3* (10) 

where tj, t 2 , and t 3 are respectively values of t calculated for 

3 3 3 ~~7 

32 , 24 , and 16 grid points at a fixed u^, and t^ is corrected value 

“7 

of t at that u^. 

In order to study the processes occurring in sheared turbulence, it is 
convenient to break the instantaneous velocities and pressures in Eqs. (1) and 
(3) into mean and fluctuating components; thus, set u. a u^ and p 

a P + p. Then, if averages are taken, and the averaged equations are sub- 
tracted from the unaveraged ones, we get, for uniformly-sheared homogeneous 
fluctuations 


8u^ dU^ 

TT " "Si ^2 


dU 2 3 u ^ 

dx^ 2 axj^ 





_ i£_ 

3X. 



3Xk 3X^^ 


( 11 ) 
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3X|^ 3X|^ 0 X 2 aXj 3X|^ 9X^ 


For Eqs. (11) and (12) the initial and boundary conditions given by Eqs. (4), 
(8), and (9) become 


u. = y a? cos 

’ if ’ 


r-t. 


.2.*b, ■ •“i’xj.b.' 

vj J J J 


<»>x„2.»b. ■ <”>, 
J J 


for any b.. 

J 

From Eqs. (11) and (12) written at a point P and similar equations writ 
ten at a point P', we can construct the following two-point correlation equa 
tions for uniformly-sheared homogeneous fluctuations^*^: 


g ouj au^ auj g 

— TylJ = -6-1 ^ ^ « jl ^ " dTJ ''2 TFJ- ^ 


— ^ (u .u'.u,' - U .U.IjT) ^ pu'. 
' 1 J k jr. ^ j 


3 u"iU^ 

TFT ^i^ * ^ 

J II 


J 1 

3rj^ ar. 


dUj 3 U.U 2 3 U^Uj^U^ 
dx^ 3T^ ar^^ ar. ’ 
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(1 


and 
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, *1 Wj 

•'■k * ^ kri ■ •''k ■ 


( 18 ) 


where the unprimed and primed quantities are measured at P and P', 

respectively, and - x^. Equations (16) to (18) show that the two 

point correlations for uniformly-sheared turbulence are functions only of t 

and not of t. That is, the turbulence is homogeneous. That is not the case 

for nonuniformly-sheared turbulence, for instance for Uj^ proportional to 
2 

instead of to x^. 

A single-point correlation equation for u.u . can be written as 

* J 


dU 


1 


— r u.u.- -6 ., -1 — U .U« 
3t 1 J il dx2 J 2 


- « 
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3X. 3X. 


(19) 


Equations (16) to (19) are useful for studying the processes occurring in uni- 
formly-sheared turbulence. 

THE LINEARIZED PROBLEM 

Equations (11) and (12) are linearized by neglecting the terms 
2 

3( u.U|^)/ 3X|^ and a (U|^u^)/3X|^ dx^^^. The numerical solution, with initial and 
periodic boundary conditions given by Eqs. (13) to (15), then proceeds as in 
the nonlinear case. 

We can obtain an analytical solution for unbounded lineariied fluctuations 
by using unbounded three-dimensional Fourier transforms as in Ref. 10. The 
solution does not satisfy constant periodic boundary conditions. Instead of 
working with the averaged equatiois used in Ref. 10, it is instructive to work 


8 




ORIGINAL PAGE IS 
OF POOR QUALITY 

with the unaveraged ones, and use the initial condition given by Eq. (13), 

instead of those in Ref. 10. In this case the Fourier transforms must be gerv- 

eralized functions (a series of delta functions), but the method of solution 

is the same as that in Ref. 10. Since U 2 is the velocity component which 

11 12 

is critical in maintaining the fluctuations against the dissipation, * we 
will, for brevity, calcu"late only that component. Equation (11) for U 2 
and Eq. (12), when linearized, are independent of u^ and u^. The solu- 
tion obtained by using the initial condition (13) is 


3 

u” cos - aqJ'tX 2 ^ 

p = / sin - aqJtX 2 ^ 

rt 1 


where 


^ 

2- V 


n n 
82 q 

2 — 

,,„n„n. . 2„n ,2 

2aq^q2t + a q^ t 


exp 




n n. . 1 2 ... 

^*^1^2^ y 3 dl 




P'’ = 


-2a 


n n n 


(2 ,2 X 

- 2 aq^q 2 t + a q^ t ^ 


exp 


-t^q” - aqjqjt + j a^qj t^^ 


(20) 


( 21 ) 


( 22 ) 


(23) 


2 2 2 2 

a = dUj/dx 2 , q" = qf + qj ^ ^iven in the 

initial conditions (Eqs. (13) and (6)). Mean values are obtained by 
integrating over all space. For instance. 


PUo = 
c 



1 

7 



(24) 
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It is clear from the form of Eqs. (20) and (21) that the solution does not 
satisfy constant periodic boundary conditions. By omitting the term 
-(dU^/dX 2 )x 2 au^/aX| as well as the nonlinear terms in Eqs. (11) and (12), we 
can, however, obtain a simple analytical solution which satisfies those 
conditions: 



This solution is useful for checking the numerical calculations and for study- 
ing the effect of the term (dUj/dx2)x23U2/ax^ on the fluctuations. 

For discussing the linearized case for constant periodic boundary condi- 
tions, it is convenient to convert Eqs. (11) and (12) to a spectral form by 
taking their three-dimensional Fourier transforms. This gives, for U 2 , on 
neglecting nonlinear terms. 



where 


n,»v 

i,(r) 





(28) 
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1 


or 
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? is the wavenumber vector, and f2 is the Fourier transform of U2* 

Note that a finite transform is used in the x2"direction in order to satisfy 
periodic boundary conditions at X2 ■ -t,x. The initial condition on 
f 2 is given by 



1 n 

7 ^2 « n 
*' 2‘’2 




( 30 ) 


where 6 (x) is the Dirac delta function 


(J ‘(‘i - ■ '(oi")) 

and 6 is again the Kronecker delta ja ^ f(<2) “ ^(^2))’ 

\x2q2 V f) 

RESULTS AND DISCUSSION 

As mentioned earlier, the U2~component of the velocity fluctuation (in 

the direction of the mean velocity gradient) is crucial in maintaining the 

11 1 ? 

turbulence against the dissipation. ’ Therefore when, for brevity, only 
one component of the velocity fluctuation is discussed, that component is 
chosen as U2. More will be said about the maintenance of tne turbulence 


later. 


1/2 


Figure 1 shows the evolution of ^2/0^ at a fixed point in space for 
a high Reynolds number, as calculateo from the full nonlinear equations. In 


11 




spite of the nonrandom Initial conditions (£qs. (13) and (6)), the velocity 

fluctuations have the appearance of those for a random turbulence. The dashed 
1/2 
~7 

curves for u^/Pq are for Initial conditions perturbed approximately 
0.1 percent. The perturbed curves at first follow the unperturbed ones but 
eventually depart sharply. Although the appearance of the curves In Figs. 1(a) 
and In 1(b) differs considerably, the perturbed curves In the two figures take 
about the same length of time to break away from the unperturbed ones. A very 
small perturbation of Initial conditions causes a large change In the values 
of U 2 except at small times. On the other hand the root-mean- square 
values of the velocities change smoothly with time and are unaffected by the 
perturbation of the inital conditions. These features are characteristic of 
turbulence. (Although the root-meari- square curve in Fig. 1(a) appears ho» i- 

zontal, it eventually goes smoothly to zero when extended.) 

1/2 

~7 

A striking feature of the curves for U 2 /'Uq in Fig. 1 is the small 
scale structure exhibited for sheared turbulence (Fig. 1(b)) when compared with 
the structure for no shear (Fig. 1(a)). This shear-related small-scale struc- 
ture is produced by the term -(cRJi/dx 2 )x 2 au^/axi in Eq. (11), from which 
is obtained the term -(dUj/dx 2 )r 2 aTTjuy 9r^ in the correlation Eq. (16). 

If we take the Fourier transform of that term we get the mean-gradient trans- 
fer term in the spectral equation corresponding to Eq. (16).^^ Its effect 
in transferring energy to small-scale components is similar to that uf the 
nonlinear transfer term (the Fourier transform of the triple-correlation term 
in Eq. (16)). The production of the small-scale structure by the shear might 
be thought of as due to a stretching of the random vortex lines in the turbu- 
lence by the mean gradient. 
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Although we first discussed a mean-gradient transfer term in 1961,^® the 
present results give the first graphic demonstration of the ef 'ectiveness of that 
term in producing a small-scale structure in turbulence, jince that is a linear 
effect, we can study it either by the full nonlinear solutions already considered 
in Fig. 1 (which contain linear as well as nonlinear effects), or by linearized 
solutions. Velocity fluctuations obtained by the latter are plotted in Fig. 2 
(Eq. (20)). The presence of small structure in the curves for dUj^/dx 2 » 4434, 
and its absence in those for dUj/dx 2 = 0 are apparent. The curve for no 
shear decays monotonically to zero when extended. This is in contrast to the 
nonlinear case in Fig. 1(a) for no shear, where at least larger fluctuations 
are present. The linearized curves for dUj^/dx 2 = 4434 in Fig. 2 follow 
closely the nonlinear ones in Fig. 1(b) for small times. Likewise the linear- 
ized curves in Fig. 2 for periodic boundary conditions follow closely those 
for unbounded conditions for small times. For larger times the fluctuations 
for unbounded conditions continue to decay, whereas those for constant peri- 
odic boundary conditions grow. The development of structure in the curves for 
unbounded conditions is produced by the term aq^tX 2 in the argument 
of the cosine in Eq. (20) (a = dU^/dx 2 ). This term arises from the term 
-ax 2 3 U 2 / 3 Xj in Eq. (11), as is evident from its absence in Eq. (25), 
where the term - ax 2 8 U 2 /axj has been neglected. For constant periodic 
boundary conditions, small-scale structure in the fluctuations or the transfer 
of energy between wave numbers is produced by the term containing the summation 
over <2 ^9* (27). That term is the Fourier transform of -ax 2 3Up/3x, 

in Eq. (11). From its form we see that it can produce a complicated inter- 
wavenumber interaction. The quantity at each interacts with 9 ^ 
at every other allowable < 2 ' Evidently a difference between the solutions 
for unbounded conditions and those for constant periodic conditions is that 


1 


f 


only fluctuations at integral values of are possible when periodic 
conditions are imposed, whereas for unbounded conditions, fluctuations are 
possible at all values of k ^. 

Although the linear term -ax^ 3U2/3Xj is effective in producing 
oscillations, even in the absence of nonlinear effects (Fig. 2), the curves 
lack the random appearance of those in Fig. 1(b). Evidently the only way we 
can have a linear turbulent solution is to put the turbulence in the initial 
conditions, as in Ref. 10. Both -ax 2 3U2/3Xj^ and the nonlinear terms 
in Eq. (11) are necessary to produce the small-scale turbulence in Fig. 1(b) 
from nonrandom initial conditions. The former acts like a chopper which chops 
the flow into small-scale components. While the latter also do that, their 
most visible effect here is to produce randomization. As in Ref. 1, the ran- 
domization might occur as a result of the presence of strange attractors in 
the flow, by proliferation of eddies or harmonic components (with the loss of 
identity of the individual eddies), or by both (see Kef. 1 for a discussion of 
these possibilities) . 

According to Eq, (20), the manufacture of small-scale fluctuations takes 
place only in the X 2 - direction. Figure 3 shows how this has taken place 
at a moderate time for a linearized case. A similar plot for the nonlinear 
case is shown in Fig. 4. The randomizing e^^fect of the nonlinear terms is 
evident. ^^2 

~T 

Figure 5 shows li^/Uq for the nonlinear case, plotted against Xj, 
rather than against X 2 as in Fig. 4. The curves show some development of 
small-scale structure in the Xj^- direction due to the interaction of the 
directional components in the nonlinear case. For the linearized flows devel- 
opment of structure occurs only in the X 2 - direction. 
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Cross-correlation coefficients u^Uj/u^ Uj (W J) are plotted 

2 

against dimensionless time for the nonlinear case In Fig. 6. Although 

~T 

a at t « 0, the Initial cross correlations are not zero but are all 

positive and equal. However because of the apparent randomization of the flow 
u^uj and l^uj approach zero as time Increases. On the other hand the values 


of the turbulent shear stress u^U 2 change from positive to negative and 

remain negative because of the dynamics of the Imposed mean shear. The 

presence of the mean velocity gradient dU]^/dx 2 causes Uj^ to be likely 

negative when Is positive, so that uju^, the correlation between the 

two, is negative. The waviness in the curves in Fig. 6, as well as that in 

some of the curves in later figures (e.g.. Fig. 9), is probably caused by non- 

random structure in the flow, possibly that produced by the linear term 

-(dUj/dx 2 )x 2 3u./8Xj in Eq. (11) (Fig. 2). 

The evolution of the mean-square components of the velocity fluctuations 

is plotted in Fig. 7. After an initial adjustment period all of the compe- 

ls 

nents increase with time, in agreement with experiment and the 

numerical results in Ref. 3. The numerical results in Ref. 4, on the 
~7 7 ~ 

other hand, show U 2 and u^ decreasing at all times. Those results are 

apparently related to the use of periodic boundary conditions on a sheared 

numerical grid rather than on a stationary one. Reference 3 also used a 

sheared grid but applied a correction (the grid was remeshed), so that the 

results are similar to the present ones in Fig. 7 for a stationary grid. 

-7 

Our u^ - component is the largest of the three, U 2 is smallest, and 

~7 ~7 13 

u^ lies slightly above u^, in agreement with experiment and previous 


numerical results. 


-7 


The effect of discretization error on the numerical results for U 2 is 

3 3 3 

shown in Fig. 8. Curves are plotted for 16 , 24 , and 32 grid 


points, together with a fourth-order extrapolation to zero grid-point spacing 
(an infinite number of grid points). The differences between the results for 

3 

32 points and the fourth-order extrapolation are small, but increase some- 
what at large times. 

Figure 9 shows the evolution of velocity pressure-gradient correlations. 
(Parts of some of the curves are omitted to avoid confusion.) The velocity 
pressure-gradient terms in the one-point correlation equation (Eq. (19)), to- 
gether with the production terms, are responsible for maintaining the turbu- 
lence against the dissipation (given by the last two terms in Eq. (19)). 

~7 ~~7 

There are no production terms in the equations for 302/ at and 3u^/3t 

/ \ ~7 ~7 

(*il *’^1^**2 *jl “i^2 zero! . Thus U2 and u^ 

generally receive energy only from the Uj^ - component, whose equation has 

a nonzero production term. Equation (19) shows that in oroer to do that, 

u. 3P/3X. and u. ap/3x. must be negative for i = j = 2, 3 and positive 

J ' ^ J 

for i = j = 1. Figure 9 shows that is actually the case for constant periodic 

boundary conditions except for an initial adjustment period, so that th<> turbu- 

lence is maintained (Fig. 7). The maintenance of the U2 - or u2-compoi ant is 

particularly critical because if U2 goes to zero, the Reynolds shear stress 


Uj^u^ in the production term of the Uj^ - equation (see Eq. 19)) will go tc 
zero and there will be nothing to keep the turbulence going. All the compo- 
nents will then eventually decay. That is what happens in the linearized 
analysis for unbounded turbulence in Fig. 9 and Ref. 10. 

A comparison between the nonlinear results for U2 and various linearized 
solutions is given in Fig. 10. The same initial conditions are used for all 
the cases (Eqs. (4) or (13) and (6)). For all of the results, except those for 
the unbounded linearized case (obtained by using unbounded Fourier transforms 
(Eq. (20)), the crucial U2 component eventually increases, so that the 


lb 


turbulence or fluctuations are maintained. In the unbounded linearized case 

~7 ~7 

u^ decreases at all tinv»s. That was expected, since the U2 - results for 

that case in Refs. 10 and 11 (for different initial conditions) decreased 

at all times. Somewhat unexpected are the linearized results for constant 

periodic boundary conditions, which show that the fluctuations a:*: maintained 

for those cases. Whereas Fig 9. shows that in the unbounded case the 

-7 

velocity pressure-gradient correlations remove energy from the U2 - component 

and cause the fluctuations to decay as in Ref. 10, the imposition of constant 

periodic boundary conditions changes the sign of those correlations and brings 
—7 

energy into U2, so that the fluctuations are maintained. Equation (25), 

which satisfies periodic boundary conditions, shows that, at least when the 

—7 

term -(dUj^/dx2)x2 au^ax^ in Eq. (11) is neglected, U2 increases at large 
n n n^ 

times if 2aq”q2 > q" for at least one n. 

Comparison of the linearized case for periodic boundary conditions in 

Fig. 10 with the corresponding nonlinear case shows that the nonlinear terms 

~7 

have a stabilizing influence. That is, the values of U2 increase more 

slowly for the nonlinear case. Moreover comparison of the curve for the 

linearized case with periodic boundary conditions and with the term 

-(dUj/dx2)x2 in Eq. (11) missing (Eq. (25)) with the corresponding 

curve for that term included shows that the presence of that term also has a 

stabilizing influence. Since neglect of that term is equivalent to neglecting 

10 

-7 

the mean-gradient-transfer term in the spectral equation for U2 , we can 

consider the latter term as stabilizing. Thus both the nonlinear transfer term 

associated with triple correlations and the linear mean-gradient transfer term 

~7 

in the spectral equation for U2 are stabilizing. The reason is that both 
terms transfer energy to small eddies where it is dissipated more easily. 
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It is of interest that the onfr-point correlation equation for au^Uj/at 
(Eq. (19)) contains neither a term associated with velocity-gradient transfer 
nor with nonlinear transfer. That is, both of those processes give zero 
direct contribution to the rate of change of u-u^; they only change the 
distribution of energy among the various spectral components or eddy sizes. 

This spectral transfer, of course, still affect the way in which u^u^ 
evolves (see Fig. 10). Even though Eq. (19) contains no transfer terms, the 
transfer of energy among the various spectral components of the velocity alters 
the terms that do appear in Eq. (19), so that du -u /at is affected indirectly. 

* J 

That is not a small effect. 

The modified linear solution given by Eqs. (25) and (26) (dash-dot-dot 
curve in Fig. 10) is the simplest solution in which the fluctuations can be 
maintained against the dissipation. In obtaining it the only mean-gradient 
term retained in the equations for U 2 (Eqs. (11) and (12), i = 2) is 
-2(dUj/dx2)au2/ axp a source term in the Poisson equation for the pressure. 

If that term is also neglected, u^ decays and, as discussed earlier, all of 
the components of the fluctuations decay. Moreover, as shown in Fig. 10 and 
already discussed, the term -(dU^/dx 2 )x 2 3u^/3x^ in Eq. (11) is stabilizing, 
so it is of no help in maintaining the fluctuations. Thus, at least in the 
linearized case, the presence of the source term -2(dU^/dx2)3U2/ 3x^ in the 
Poisson equation for the pressure is necessary for maintaining the fluctua- 
tions. That term should play a similar important role in the maintenance of 
nonlinear turbulence, although in that case it is hard to separate the linear 
effects from the nonlinear ones. In particular, the role of the nonlinear 
source term in the Poission equation for the pressure remains unclear, al- 
though it may have an effect similar to that of the linear source term. 
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Figures 11 and 12 show the approach to Isotropy of nonlinear uniformly 


sheared turbulence when the shear is suddenly removed. Although the shear 

-7 

produces considerable anisotropy^ the components u^ of the mean-square 

fluctuation approach equality upon removal of the shear and remain equal. The 

velocity pressure-gradient correlations in Eq. (19) are thus successful in 

transferring energy among the various directional-components in such a way that 
-7 -7 

eouality of the u^ is produced. We note that u^ continues to increase 

for a short time after the shear is removed, probably because it receives 

-7 -7 

energy from both Uj and u^. 


In addition to equality of the u. , zero cross-correlations u.u. are 

1 1 0 

required for isotropy. Figure 12 shows that u^u^, which is nonzero when the 
turbulence is sheared, approaches zero when the shear is removed, and along 
with the other cross-correlations, remains close to zero. The destruction 


of Uj^U 2 > apparently by nonlinear randomization effects, occurs over a finite 
time period rather than instantaneously on removal of the shear. 

Another expected effect of removal of the mean shear is that the small 
scale structure produced by the term -(dUj^/dx 2 )x 2 au^/ax^ in Eq. (11) should 
die out. According to Fig. 13, that occurs almost immediately when dUj^/cx 2 
goes to zero, evidently because of the large fluctuating shear stresses between 
the small-scale eddies. 

CONCLUDING REMARKS 


According to our results for both sheared and unsheared flow at higher 
Reynolds number, the nonlinear structure of the Navier-Stokes equations is 
such that an apparently random turbulence can develop from nonrandom initial 
conditions. The presence of a mean gradient, in addition to producing a non- 
zero turbulent shear stress, produces small-scale fluctuations in the flow. 
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These can be attributed to a mean-gradient transfer term in the spectral equa- 
tion for the velocity fluctuations,^^ or to an equivalent term in Eq. (16) 
or (11). However the small-scale fluctuations produced by that term alone 
(linear solution) are essentially nonrandom. Evidently the only way we can 
have a linear turbulent solution is to put the turbulence in tne initial con- 
ditions, as in Ref. 10. In order to produce the observed small-scale turbu- 
lence from nonrandom initial conditions, the presence of both the linear mean- 
qradient transfer term and the nonlinear terms in the equations is necessary. 
The former term, or its equivalent in Eq. (16) or (11), acts like a chopper 
which chops the flow into small-scale components, and the latter terms, while 
tliey also produce small-scale components, act most visibly here as randomizers. 

In all of the cases calculated with constant periodic boundary conditions, 
including both linear and nonlinear flows, the velocity pressure-gradient cor- 
relations are successful in distributing the energy among the directional com- 
ponents (Fig. y), so that tne turbulence or the fluctuations are maintained. 
This is in spite of the presence of a production term in the equation for only 
one of the components. Botn the linear mean-gradient transfer term and the 
nonlinear terms mentioned in the last paragraph have a stabilizing effect. 

That is, they cause the fluctuations to increase at a slower rate. It is 
snown that, at least for the linearized solution with constant periodic boun- 
dary conditions, a meafv-gradient source term in the Posisson equation for the 
pressure is necessary for maintaining the fluctuations against the dissipa- 
tion. That term snould play a similar important role in ttie maintenance of 
nonlinear turoulence, altliough in that case it is hard to separate the linear 
effects from the nonlinear ones. In particular, the role of the nonlinear 
siH.rce term in the Poisson equation for the pressure remains unclear, although 
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it may have an effect similar to that of the linear source term. For the lin- 
earized unbounded solution (obtained by using unbounded Fourier transforms) 
the fluctuations decay, as expected from the results of Ref. 10. 

When the mean-velocity gradient is suddenly removed, the turbulent shear 
stress goes to zero in a finite time period, and the velocity pressure- 
gradient correlations cause the turbulence to attain the isotropic state. The 
intensities of the directional components become and remain equal. In addi- 
tion, the small-scale structure produced by the mean-gradient transfer term 
quickly vanishes. 
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